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Abstract 

A random walk model in a one dimensional disordered medium with an os- 
cillatory input current is presented as a generic model of boundary perturba- 
tion methods to investigate properties of a transport process in a disordered 
medium. It is rigorously shown that an admittance which is equal to the 
Fourier-Laplace transform of the first-passage time distribution is non-self- 
averaging when the disorder is strong. The low frequency behavior of the 
disorder-averaged admittance, (x) — 1 ~ where \i < 1, does not coincide 
with the low frequency behavior of the admittance for any sample, x ~ 1 ~ UJ - 
It implies that the Cole-Cole plot of (x) appears at a different position from 
the Cole-Cole plots of x of any sample. These results are confirmed by Monte- 
Carlo simulations. 

PACS numbers: 02.50.Ey, 05.40.Fb, 72.80.Ng 
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Response of a system to an external field is a standard information to be utilized in 
the study of condensed matter. Recently, the frequency response method and the in- 
tensity modulated photocurrent spectroscopy 0| have been introduced, where an oscillatory 
perturbation is applied at one end of a system and a response from the other end of the 
system is measured. Thus, these experimental techniques can be considered to belong to 
a generic method which can be called the boundary perturbation method (BPM). In the 
presence of a periodically forced boundary condition on one end of a system, the output 
from the other end of the system is in proportion to the perturbation in the linear regime 
and the proportionality constant is called the admittance. The frequency dependence of 
the admittance contains various information of the dynamics of the system. The BPM is 
expected to provide useful information on transport properties of disordered media. 

The analysis of transport of particles in a random medium attracts great interest since 
the transport mechanism is basic for the understanding of many physical phenomena, from 
electrical conductivity to thermal properties ||. Usually it is assumed that a sample used in 
experiments is sufficiently large that it can be considered to be composed of a large number 
of sub-systems. If each of sub-systems is itself macroscopic, the boundary effect is negligible 
and each sub-system may be considered to be a realization of the system with a particular 
choice for the disorder. Thus, a measurement of any observable in such a system corresponds 
to an average over all the sub-systems, i.e. an average over the ensemble of all realizations 
of the disorder. Systems for which this assumption is valid are said to be self-averaging 
(SA). Since the assumption implies no sample dependence, it guarantees reproducibility of 
experimental results in any sample. 

However, it is known that the disorder-averaged quantity follows a dynamics that is 
different from the dynamics of the quantity of each sample. If the dynamics of particles 
in a disordered sample is assumed to be stochastic and Markovian, i.e. to follow a master 
equation, the disorder-averaged master equation has a memory effect |4j]. Thus, the dynamics 
of an averaged quantity may be non-Markovian even though the process for each sample is 
Markovian. There are several reports on absence of SA in transport phenomena in disordered 
media: for mean square displacement ||, for mean first-passage time of the Sinai model ||. 

In this letter, we study a random walk model in a one dimensional disordered medium 
with an oscillatory input current as a generic model of the BPM. It is rigorously shown that 
the admittance, which is equal to the Fourier-Laplace transform of the first-passage time 
distribution (FPTD), is non-self-averaging (non-SA) when the disorder is strong. 

We consider a general random walk in a one dimensional lattice segment of N + 1 sites. 
The lattice sites are denoted by integers, n = 0, 1, . . . , N. The dynamics of a random walking 
particle can be described by a master equation for the probability P n {t) that the particle is 
at site n at time t > 0. The master equation for the model is written as 

P n (t) = W n n -iP n -i(t) - (W n -l n + W n+ i n )Pn{t) + Wn n+1 P n+1 (t), (0.1) 

where w m n denotes the random jump rate of a particle from site n to site m. The probability 
distribution of jump rates characterizes the random medium. We introduce a perturbation 
at the left end, so that the equation for Po(t) is given by 

Po(t) = -wi oP Q (t) + w ^(t) + J(t). (0.2) 
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J(t) is the oscillatory current perturbation at site 0. The right end of the system is assumed 
to be a sink and the equation for Pjy-i is given by 

PN-l{t) = W N -i N~2PN-2{t) ~ {w N -2 N-l + ™N iV-l)^V-l(*)- (0-3) 

Since the current perturbation J(t) oscillates in time around a positive average with the 
amplitude A J, the response of the output current wn jv-i-P/v-i(£) oscillates around its 
stationary state with the amplitude wn jv-iAPjv-i at the same frequency with a phase- 
shift. The admittance is defined as the ratio of these two amplitudes 

X{w) = • (0.4) 

We first note that the admittance can be related to the first passage time distribution 
Fn o(t) (FPTD), which is the probability density that the particle which starts at site at 
time arrives for the first time at site N at time t. Since site N is a sink, the FPTD is given 
by the output current from site N when there is no input current and the particle starts 
from site at time 0. The Fourier-Laplace transform of a system of master equations for the 
case is the same as a system of equations for AP n /Aj derived from the master equations 
Eqs. ( p.l| . p.2|JX3l) . It implies that F N (s — iuS) for the case where there is no input current 
and the particle starts from site at time is identical to wn jv-iAP/v-i/A J. Thus, the 
admittance is equal to the Fourier-Laplace transform of the output current in the problem 
described above, i.e. the FPTD. 

We make the low-frequency expansion of the admittance to see the behavior near the 
static limit. Since the admittance is given by the Fourier-Laplace transform of the FPTD, 
the admittance at zero frequency equals to one due to normalization of the FPTD. The 
mean first-passage time (MFPT) i which is the first moment of the FPTD is given by 
i — i dx(uj = 0)/cLj. Thus, the low-frequency expansion of the admittance is given by 

X {uj) = 1 -i t uj + 0(uj 2 ). (0.5) 

Since the MFPT is given by [7J 

N-l 1 N-2 x N-l i . 

*'=£— —H (o-6) 

fc=0 W k-l k k=0 W k+l k i=k+1 j=k+1 W j+l j 

the MFPT t is larger than J2k=o V^fc-i k- By taking the average of Eq. (p.6|) over the 
distribution for the jump rates, one sees that the average of the MFPT (t) is larger than the 
first inverse moment of a random jump rate (1/w). Thus, when the first inverse moment of 
a jump rate diverges, which is called strong disorder, the disorder-averaged MFPT diverges. 
It indicates that the disorder-averaged admittance is non-analytic. Since the admittance is 
a function of iu, we can write as 



(xH) - 1 ~ (iu;)", (0.7) 

where < /i < 1. Since this behavior of the disorder-averaged admittance is completely dif- 
ferent from the behavior of the admittance of each sample given by Eq. ( p.5|) , the admittance 
is non-SA when the disorder is strong. 
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FIG. 1. Low- frequency behavior of the imaginary part of the admittance when the exponent 
of the probability distribution of jump rates is 1/2 and the length of a chain is 100. Each solid 
line represents the admittance for each of 50 samples. The dashed line represents the admittance 
averaged over 5000 samples. Though the admittance of each sample is proportional to u>, the 
averaged admittance is proportional to \fuj. 



On the other hand, when the first inverse moment of a jump rate exists, it is easily 
shown from Eq. ( p.6|) , that the MFPT is finite provided that there is no correlation between 
Wk+i k and w i+ i i; where i > k + 1. The condition of independence of the neighboring jump 
rates holds for the site disordered model called random trap model and the bond disordered 
model called random barrier model. Thus, when the first inverse moment of a jump rate 
exists, the low frequency behavior of the disorder-averaged admittance is as same as the one 
of each sample given by x ~ 1 ~ Furthermore, if the second inverse moment of a jump 
rate exists, the prefactor of u does not depend on a realization of the chain when the chain 
is sufficiently long. 

In order to test the foregoing observation, we performed a computer simulation in one 
dimensional chain of 100 sites. Figure [l] shows the low frequency behavior of the imaginary 
part of the admittance obtained by the Monte-Carlo simulation where the relation Wi + i j = 
Wi i + \ is assumed and the probability distribution of jump rates is given by 

P( w ) = ! anja ~ l if < w < 1 (0 8) 

1 otherwise, 

where < a < 1. The power-law probability distribution corresponds to the waiting-time 
distribution ^(t) ~ used in the study of dispersive transport where a is proportional 

to temperature ||, and it has been shown that the power-law distribution is common for 
activation processes with random activation energy ||. 

In the literature of experiments of the BPM [0,0, the Cole-Cole plot of the admittance 
is employed for the analysis. It is a parametric plot of the imaginary part of the admittance 
against the real part. The shape of the Cole-Cole plot is independent of the length of the 
chain for a sufficiently long chain due to time-scale invariance of the Cole-Cole plot. Thus, 
the Cole-Cole plot shows the size-independent characteristics of a medium. 

In order to see the non-SA property of the admittance in the Cole-Cole plot, we first 
note that the inverse of the admittance satisfies the following recursive relation: 
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-1/ \ . IUJ w n—l n\ -if \ w n-l n —1 / \ /n n\ 

Xn M = ( + 1 + )X»-l( W ) X„- 2 M, (°- 9 ) 

n W n +i n Wn+1 n 

where x~ 1 (c<j) is the inverse of the admittance for a chain of length n [10]. The recur- 



sive relation proves inductively that the inverse of the admittance obeys the following four 
inequalities in the frequency region [0, oj n ], where uj n is the positive smallest zero of Xn 1 ''- 

— 1 " n 

Xn > 

Xn _^ >0 _x n (0.10) 



< 



Xn ^ Xn—1 
-1' -1 ' 

, Xn ^ Xn—1 • 



Since the second inequality shows u n < u n -\, the fourth inequality implies x n V < Xn-i < 
• • • < x^f 1 = 1 in the region [0, uo n ]. It is straightforward to show that x 1 ' < 1 implies that 
the Cole-Cole plot of the admittance of each sample can exist only outside a semicircle whose 
center is at (1/2, 0) and radius is 1/2, which we call the Debye semicircle, when uo < uj n . 

On the other hand, we can show in the following way that the Cole-Cole plot of the 
averaged admittance is located inside the Debye semicircle when the disorder is strong. We 
consider the angle 9 between the tangent of the Cole-Cole plot at (1, 0) and the horizontal 
axis. For Eq. (|0.7p , 9 = /U7r/2, which is less than n/2 since \i < 1, is obtained. However the 
angle 9 for the Debye semicircle is equal to tt/2. It implies that the Cole-Cole plot of the 
averaged admittance appears inside the Debye semicircle when the disorder is strong. Since 
the Cole-Cole plot of the admittance of any sample appears outside the Debye semicircle, 
the Cole-Cole plot clearly shows non-SA property of the admittance. Figure ^| shows the 
Cole-Cole plot of x an d (x) obtained by Monte Carlo simulation when a = 1/2 in Eq. 
(p.8[). Since the admittance of each sample at the same frequency scatters outside the Debye 
semicircle, the admittance averaged at the same frequency is inside the Debye semicircle. 

In the case that the first inverse moment of a jump rate exists, the angle 9 of the Cole- 
Cole plot of the averaged admittance is tt/2, which is as same as the Cole-Cole plot of each 
sample, because of the regular behavior as (x) — 1 ~ iuj. 

In conclusion, we have presented a generic stochastic dynamical model in one dimensional 
medium of the BPM and have shown that the admittance has important information about 
transport in disordered media, e.g. the FPTD. Thus the BPM will be a powerful technique to 
investigate the dynamical process in random media. We have also shown that the admittance 
is non-SA when the disorder is strong. Since the absence of SA is due to non-analyticity 
introduced by an average over infinite number of samples, the behavior of the averaged 
admittance over infinite number of samples given by a theoretical analysis does not coincide 
with that of the admittance observed in experiments. 

Non-SA behavior, i.e. strong sample to sample dependence may be observed in other 
experiments. For example, it is known that anomalous system-size dependence of the mobil- 
ity, /i ~ jV 1-1//Q , is observed in the time-of-flight experiment of amorphous semiconductors 
H. The anomalous system-size dependence is due to the fact that the smallest value of 



random jump rates depends on the length of a chain Since the smallest value domi- 
nates the mobility, the mobility is also a random variable. Thus, the anomalous system-size 
dependence implies absence of SA. 
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FIG. 2. The Cole-Cole plot of the admittance when the exponent of the probability distribution 
of jump rates is 1/2 and N = 100. Each solid line represents the admittance for each of 30 samples. 
The dashed line represents the admittance averaged over 5000 samples. The circles represents the 
admittance of each sample at to = 5 x 10 -7 . Although the circles scatter outside the Debye 
semicircle, the diamond, i.e. the averaged admittance is inside of Debye semicircle. 

It is important to note that the system treated in the present letter is purely one- 
dimensional and it is still an open problem whether the admittance for higher dimensional 
media are SA. Since one can consider a diffusing particle's path a one dimensional chain and 
a higher dimensional medium can be regarded as a bundle of many realizations of a particle's 
path, the total output from the medium is given by the sum of one dimensional currents. 
Thus, the admittance for such a system may be SA, since the number of paths is very large 
when the medium is macroscopic. Thus, in order to see the non-SA behavior described here 
one needs to work in one dimensional system or the case where the translational invariance 
is violated in only one direction since there is still possibility that the SA assumption is valid 
for two or three dimensional systems. 
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tion, Science, Sports and Culture. One of us (M.K.) is thankful to the Institut fur 
Festkorperforschung, Forschungszentrum Jiilich for its hospitality, where part of the present 
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